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'Gravitational memory' refers to the possibility that, in cosmologies with a time-varying gravi- 
tational 'constant', objects such as black holes may retain a memory of conditions at the time of 
their birth. We consider this phenomenon in a different physical scenario, where the objects under 
consideration are boson stars. We construct boson star solutions in scalar-tensor gravity theories, 
and consider their dependence on the asymptotic value of the gravitational strength. We then 
discuss several possible physical interpretations, including the concept of pure gravitational stellar 
evolution. 
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I. INTRODUCTION 

A few years ago, Barrow introduced the concept of 
gravitational memory pj by posing the following prob- 
lem: what happens to black holes during the subsequent 
evolution of the universe if the gravitational coupling G 
evolves with time? He envisaged two possible scenarios. 
In the first, dubbed scenario A, the black hole evolves 
quasi-statically, in order to adjust its size with the chang- 
ing G. If true, this means that there are no static black 
holes, even classically, during any period in which G 
changes. In the alternative possibility, scenario B, the 
local value of G within the black hole is preserved, while 
the asymptotic value evolves with a cosmological rate. 
This would mean that the black hole keeps a memory of 
the strength of gravity at the moment of its formation. 
Further analysis of the striking phenomena which arise 
in both of these scenarios was made in Ref. Black 
holes are of particular interest in this context, because 
primordial black holes may have formed in the very early 
stages of the Universe. At those times no direct evidence 
of the strength of gravity is available; nucleosynthesis is 
the earliest epoch at which significant constraints apply. 
Gravitational memory may provide a unique probe of 
these early epochs. However, in this paper our aim is 
to consider the gravitational memory phenomenon in a 
completely stellar arena. The motivation is twofold. On 
one side, since there are no singularities or event horizons 
as in the black hole case, there are many calculational ad- 
vantages which may provide a more direct route towards 
shedding light on which of the scenarios stated above will 
happen in practice. On the other, a comprehensive un- 
derstanding of the effects of a varying-G cosmology upon 
astrophysical objects is far from complete. 

To make a proper study of a varying gravitational 
strength, it is imperative to operate within a self- 
consistent framework rather than simply writing in a 
variation 'by hand'. The most useful such framework 



is scalar-tensor theories p|, which have an action 
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Here g^ v is the metric, R is the scalar curvature, (j) is 
the Brans-Dicke field, C m is the Lagrangian of the mat- 
ter content of the system, and uj{4>) gives the strength 
of the coupling between the scalar (f> and the metric. 
The special case of a constant w is the Jordan-Brans- 
Dicke (JBD) theory, and general relativity is regained 
in the limit of large to. The gravitational constant G 
of the Einstein-Hilbert action is replaced by a dynam- 
ical field, In the usual applications, it is cither 
spatially-constant but time-varying (cosmological solu- 
tions) or spatially-varying but time-independent (astro- 
physical solutions). Our situation will be unusual in in- 
cluding both types of variation. 

Our choice of stellar object is motivated by two con- 
cerns. Firstly, it should be simple enough as to allow us 
to isolate the effects caused by a varying G. Secondly, in 
order to show the changes due to the particular gravita- 
tional theory, it should be fully relativistic, in the sense 
that its equilibrium configuration must be obtained from 
Einstein-like field equations. There is one type of stel- 
lar object which meets these criteria, though so far it is 
known to exist only as a theoretical construct. It is the 
boson star, the analogue of a neutron star formed when 
a large collection of bosonic particles (normally taken to 
be scalars) becomes gravitationally bound. Such con- 
figurations were introduced by Kaup (J] and by Ruffini 
and Bonazzola |5| in the nineteen sixties, but current in- 
terest was sparked by Colpi et al. Q who showed that, 
provided the scalar field has a self-interaction, the bo- 
son star masses could be of the same order of magnitude 
as, or even much greater than, the Chandrasekhar mass. 
This led to the study of many properties, which are sum- 
marized in two reviews |7j . It seems possible for them to 
form through gravitational collapse of a scalar field [||, 
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though little is known as yet. The observational status of 
boson stars was analyzed very recently in Ref . || , where 
it was asked whether radiating baryonic matter moving 
around a boson star could be converted into an observa- 
tional signal. Unfortunately, any direct detection looks a 
long way off. 

As far as we are aware, only three papers have stud- 
ied boson stars in scalar-tensor theories. Gunderson 
and Jensen [jlO| studied the JBD scenario, concluding 
that typically the mass of equilibrium solutions would 
be reduced by a few percent. This work was general- 
ized by Torres jnj, who looked at several coupling func- 
tions lu(4>), chosen to be compatible with known weak- 
field limit and nucleosynthesis constraints. Similar sce- 
narios were further studied in Ref. [^2]| . In this paper we 
shall study boson stars formed at different times of cos- 
mic history and in addition consider the possible phys- 
ical evolution of such objects. One might further hope 
that results found for boson stars might be representa- 
tive of other stellar candidates, especially neutron stars 
but perhaps even indicative of what might happen with 
long-lived hydrogen-burning stars such as our own Sun. 



ip(r, t) = x(r) exp [— izut] . 
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To write the equations of structure of the star, we use a 
rescaled radial coordinate, given by 

x = mr . (7) 

From now on, a prime will denote a derivative with re- 
spect to the variable x. We also define dimensionless 
quantities by 



pi 



where mpi = G is the present Planck mass. Our 
observed gravitational coupling implies $ = 1. In order 
to consider the total amount of mass of the star within 
a radius x we change the function A in the metric to its 
Schwarzschild form, 
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II. EQUILIBRIUM CONFIGURATIONS 

We begin with a review of the formalism, which can 
be found in more detail in Ref. The material from 

which the boson star is made is a complex, massive, 
self-interacting scalar field ipi which is unrelated to the 
Brans-Dicke scalar already described. Its Lagrangian is 
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The U(l) symmetry leads to conservation of boson num- 
ber. Varying the action with respect to g^ u and <fi we 
obtain the field equations: 
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The issue of the definition of mass in JBD theory is quite 
a subtle one p3| . The above gives the Schwarzschild 
mass; there are arguments that the tensor mass is more 
appropriate Jl3| , but we have verified numerically that 
for the large couplings we use the difference is negligible. 

Note that a factor $(oo) appears in Eq. ([)[). This is 
crucial to obtain the correct value of the mass, which is 
given by 



M star = M(oo) $(oo) -™ 
m 



(10) 



for a given value of m. The <E>(oo) factor allows for the 
asymptotic gravitational coupling to be different from 
that presently observed]^] With all these definitions, the 
non-trivial equations of structure are 
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where T^ v is the energy-momentum tensor for the matter 
fields (Eq. (5) of Ref. |ll|) and T, its trace. Commas 
and semicolons are derivatives and covariant derivatives, 
respectively. 

At first, we seek static equilibrium solutions. The met- 
ric of a spherically-symmetric system can be written as 



ds 2 = -B(r)dt 2 + A(r)dr 2 
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We also demand a spherically-symmetric form for the i/j 
field. Crucially, the most general ansatz consistent with 
the static metric permits a time-dependent ip of the form 
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*Note that the $(oo) factor was mistakenly forgotten in 
Ref. |0. This affects Table III of that work, where M(oo) is 
actually M(oo)/$(oo). The $(0) values and the conclusions 
extracted from it are unaffected. 
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To solve these equations numerically, we use a fourth- 
order Runge-Kutta method, for which details may be 
found in Ref. §L% . Here we shall study two kinds of 
theories. 

1. JBD theory with u> — 400. This is comparable to 
current observational limits jl4|,[l5| ■ 

2. A scalar-tensor theory with a coupling function of 
the form 

2cj + 3 = 2Bi|1-$|- Q! , 

where we choose a = 0.5 and B\ = 5. The cos- 
mological setting of this model has been studied by 
Barrow and Parsons |l(| . 

As explained in Ref. JllJ], the results for these couplings 
may be thought of as the general behavior of any scalar- 
tensor theory, by suitable expanding a general coupling 
in a Taylor or Laurent series. The positive exponent of 
this theory must be obtained from power-law couplings, 
but the behavior of these was found to be similar to the 
pure Brans-Dicke ones. 

In addition to the freedom to choose the fundamen- 
tal parameters, there are two free boundary conditions. 
One is the value of the boson field at the centre of the 
star, cr(0), which we call the central density. The sec- 
ond is the asymptotic value of the Brans-Dicke field, 
which determines the asymptotic strength of gravity. The 
other boundary conditions are fixed by demanding non- 
singularity and finite mass @jn[l. This still leaves an 
infinite discrete set of solutions with different vj, cor- 
responding to a different number of nodes in <j(x). We 
choose the nodeless solution, which is the only stable one. 
Higher node solutions are generated in Ref. Il2| . 



III. QUASI-STATIC EVOLUTION 

We first assess the likely cosmological variation of </>, 
concentrating on JBD theory. During radiation domi- 
nation, the attractor behavior is actually exactly that 
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FIG. 1. Boson star masses as a function of $(oo); squares 
are for the JBD theory and circles the scalar-tensor theory. 
The upper panel shows models with A = 100 and a(0) = 0.06. 
The lower one has A = and cr(0) = 0.1. 



of general relativity, namely a constant <fi and a cx i 1 / 2 . 
This changes with the onset of matter domination, when 
the attractor solution becomes [l7|JT^ | 



a(t) cx i( 2 ~")/ 3 ; G(t) cx r 
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where n = 2/(4 + 3oj), so G exhibits a slow decrease. 
Assuming that matter-radiation equality took place near 
the general relativity value, at z cq = 24 000f2o/i 2 , then, 
for critical density and h = 0.5, the fractional change in 
G since equality is 



G(t ) 
G(t cq ) 



= 6000~ 1/(1j 
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For oj = 400, the ratio is 0.98, so G will have changed 
value by about two percent since equality. 

Informed by the likely range of variation, we compute 
static configurations with different values of the asymp- 
totic gravitational strength, denoted <&(cjo) = 1/G. In 
Fig. 1, we see that if A and the central density er(0) 
are kept fixed, then the mass changes with a different 
$(oo). In the JBD case, the mass is a growing function 
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FIG. 2. This shows the boson mass as a function of parti- 
cle number in the two models, with squares for the JBD theory 
and triangles for the scalar-tensor one. Both have A = and 
<r(0) = 0.1. The leftmost point in each case corresponds to 
$(oo) = 0.9 and the rightmost to $(oo) = 1.1. As M(N) is 
single valued, there can be no change of stability. 

of 4>(oo), while in the more complicated second model 
there is actually a small peak in the mass at our present 
gravitational strength. That peak is due to the special 
form of the coupling function, which selects out our own 
observed coupling strength. 

Let us now suppose that evolution does drive the sys- 
tem through a series of quasi-static states, and naively 
estimate the sorts of energies involved. Fig. 1 shows that 
variation of $(oo) by a few percent can change the masses 
of configurations by a few percent. As it happens, this is 
very similar to, or even greater than, the fraction of the 
mass-energy which can be liberated by nuclear reactions 
over the lifetime of a star, a process which for most stars 
occurs on similar timescales to the cosmological ones of 
the time- varying G. 

Before analyzing this in further detail, we shall first 
show the computation for the number density of the con- 
figurations, defined by 

^-N = Q [°° a\[^x 2 dx, (17) 
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We found that, for both models under consideration, the 
number of bosons is a growing function of $(oo). Fig. 2 
shows the bifurcation diagram, number density against 
mass. There is no cusp structure in the sense of catas- 
trophe theory, which implies that there is no change in 
the stability criterion for values of G close to the present 
one @]. In fact, we found this result for a wider range 
of $(oo) than shown. 

IV. SCENARIOS 

Now we discuss possible interpretations of the com- 
putational results described above, in the same context 



which led Barrow to propose the gravitational memory 
hypothesis. We describe the possible subsequent histo- 
ries a boson star might have, after forming at a time tf 
with a mass M(tf). 

A. The gravitational memory hypothesis 

The star remains completely static, without change in 
the values of its mass or central density. Such a situa- 
tion would be reminiscent say of a virialized galaxy after 
gravitational collapse, which has become decoupled from 
the cosmological expansion of the Universe around it and 
in particular does not participate in further expansion. 
In such a scenario the mass M(tf) is a function of the 
formation time, and hence of $(oo , tf), as well as the 
central density; the star keeps memory of the value of G 
at formation. This situation clearly cannot be precisely 
correct since the asymptotic gravitational constant does 
evolve, but it may well be that the effect on the star is 
much slower even than the cosmological evolution of G, 
so that in practice the star can be viewed as static. 

In this scenario, there is the interesting feature that 
stars of the same mass may differ in other physical prop- 
erties (their radius, for example) depending on the for- 
mation time. 



B. Quasi-static evolution hypothesis 

The opposite regime has the adjustment time for the 
star being much shorter than the asymptotic evolution of 
G. If true, the star should quasi-statically evolve, either 
changing its mass or its central density or both. In this 
scheme, when intervals of time are short enough com- 
pared with the scale of cosmological evolution, the star 
can be taken as static and the solutions computed hold 
in this limit. 

An interesting subcase is pure gravitational evolution. 
This assumes that the mass of the star is preserved all 
along the quasi-static evolution, through the central den- 
sity evolving in the appropriate way. This leads to the 
remarkable conclusion that stars may be able to evolve 
even without absorption or emission of energy, and with- 
out nuclear burning. Such quiescent evolution would be 
entirely gravitational in its nature. 

Purely gravitational evolution would doubtless have 
significant consequences for stellar evolution, and indeed 
may well be quite strongly constrained. Yet within the 
quasi-static evolution hypothesis, it is actually rather 
conservative, because the other possibility is that mass 
is lost during the evolution, if dissipative processes are 
required to keep tracking the evolutionary sequence. We 
have seen that the variation in G can easily lead to 
changes in the mass of up to a few percent, which can 
be of the same order as the energy liberated in a conven- 
tional star by nuclear burning. 
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The reduction in central density seems reasonable 
when we recall the force balance in polytropic stars. 
Since gravity is reducing in strength as time goes by, the 
equilibrium configurations can become more diffuse and 
hence drop in central density. In an extreme situation, 
one might wonder whether the continuation of this pro- 
cess can in fact lead to the complete destruction of the 
star, though this may be prevented by the initial nega- 
tive binding energy of a stable boson star, since boson 
number is conserved. 

The quasi-static evolution hypothesis has an important 
difference from the gravitational memory hypothesis; in 
it, stars of the same mass are identical in all their other 
physical properties too. 

C. Feedback on the asymptotic gravitational 
constant? 



step would be numerical simulation including dynamical 
evolution, starting from a static solution and varying the 
asymptotic gravitational coupling. This looks a promis- 
ing avenue for determining which of the proposed scenar- 
ios is the correct one, and that knowledge may allow a 
test of general relativity not just at the present epoch, 
but in the distant past, through astrophysical systems. 
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Finally, one can ask whether the formation of boson 
stars can significantly influence the 'asymptotic' gravi- 
tational constant; when the stars form gravity becomes 
weaker in their interior as the <p value increases. Follow- 
ing the results of Ref . jll[] (Table III) , this change can be 
about 1% between the internal and the external G value. 
In a static configuration, the radius at which cj> finally 
approaches its asymptotic value is quite a bit larger than 
the region in which the t/j field is localized. If a very 
high density of boson stars formed, might they be able 
to reduce the gravitational interaction strength in quite a 
significant region around themselves? It is an interesting 
possibility, though unlikely in practice as the density of 
material available to make the stars is so small; boson 
stars would be expected to be separated by similar dis- 
tances to conventional stars (which could also contribute 
to the effect). 



V. DISCUSSION 

In this paper, we have not directly tried to answer the 
question as to whether or not the gravitational memory 
phenomenon exists. Rather, we have introduced a new 
framework in which it can be studied, that of stellar sys- 
tems. As a spinoff, it allows us to introduce the concept 
of pure gravitational stellar evolution. As well as their 
possible physical relevance, boson stars have the consid- 
erable advantage of being much simpler than their black 
hole analogues. They are also much simpler than neu- 
tron stars, as they are based directly on a field theory 
description. In that light, we have studied static config- 
urations in two different scalar-tensor theories, in par- 
ticular emphasizing the dependence of the mass on the 
asymptotic value of the gravitational coupling. We have 
also been able to make some preliminary statements on 
the dynamical stability of the configurations, an issue we 
leave for further study in a future publication. The next 
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